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Abstract 

We calculate the energy and hydrostatic pressure densities of a hot quark-gluon 
plasma in thermal equilibrium through diagrammatic analyses of the statistical av- 
erage, (Qfiv), of the energy-momentum-tensor operator 0^. To leading order at 
high temperature, the energy density of the long wave length modes is consistently 
extracted by applying the hard-thermal-loop resummation scheme to the operator- 
inserted no-leg thermal amplitudes (0^)- We find that, for the long wave length 
gluons, the energy density, being positive, is tremendously enhanced as compared 
to the noninteracting case, while, for the quarks, no noticeable deviation from the 
noninteracting case is found. 
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1 Introduction 



For the long wave length { A < O^gT)^ 1 ) } or soft modes in a hot quark-gluon plasma, 
the contributions of "hard thermal loops" to their amplitudes are as important as tree 
amplitudes |l|, |, [| f|. Numerous applications have been made || of the hard thermal 
loop (HTL) resummation scheme EL 131, M, the scheme which enables consistent 
evaluations of any thermal amplitudes to leading order in the coupling constant g. 

Due to the HTL resummations, the energy and pressure densities of the soft 
modes in the quark-gluon plasma deviates from the form in the noninteracting case. 
Weldon was the first who treat this issue. Subsequently, several work [0, |8], 01 have 
been devoted to the issue. The strategy for extracting the energy-momentum-tensor 
operator for soft modes in || 0, ^§ and the energy functional in || is to use the effective 



action [HI 11 1, in which all the contributions coming from HTL resummations are 



taken in. A brief account is given in || of the energy-momentum tensors derived in 

B EE- 

The purpose of the present paper is to evaluate the energy density £{j>) and 
the pressure density V(p) of the soft modes in the hot quark-gluon plasma from 
purely "diagrammatic point of view" . By "the soft mode" we mean the mode whose 
momentum p is in the range 

g 2 T«p«T. (1) 

For analyzing the hyper-soft modes p < 0(g 2 T), a treatment beyond the HTL resum- 
mation scheme is required. Our procedure is as follows: 1) We start with the energy- 
momentum-tensor operator, M1/ , which are extracted from the QCD Lagrangian. 2) 
We evaluate the energy and pressure densities through calculating the statistical av- 
erages (Ooo) and (On), respectively. 3) In evaluating these densities for soft modes, 
to leading order at high temperature, we apply the HTL resummation scheme to the 
operator inserted no-leg thermal amplitudes, (Ooo) and (On). Through the HTL 
resummations, in addition to (formally) the lowest-order diagrams, formally 0(g 2 ) 
diagrams yield dominant contributions. 

It is worth mentioning that, in the approach || [7], [|], H|, the HTL resummations 
are performed first to construct the effective action and then the energy-momentum- 
tensor operator (for soft modes) is extracted. On the contrary, our approach starts 
with the energy-momentum-tensor operator extracted from the Lagrangian and 
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then (Goo) and (On) are evaluated on the basis of the HTL resummation scheme. It 
turns out that the resultant (0^) corresponds to the thermal average of the energy- 
momentum-tensor operator found in 0j. Relationships between various expres- 
sions for T^ u , found in || [7], §, |!|, are not immediately obvious @. This is also the 
case for the relations of our (O^) to T MJ ,'s in || ]8|, 

The plan of the paper is as follows. In Sec. II, some preliminaries, which include 
the expression for extracted from the Lagrangian, are given. The contribution to 
(©/^) from formally lowest-order diagrams are evaluated in Sec. III. In Sec. IV, on the 
basis of the HTL-resummation scheme, we calculate the contributions to (O^) from 
formally 0(g 2 ) diagrams. In Sec. V, the energy of the soft modes is calculated and, 
in Sec. VI, the hydrostatic pressure is calculated. Sec. VII is devoted to summary 
and discussion. 



2 Preliminary 

We consider an SU(N) gauge theory with Nf flavors of massless quarks in the fun- 
damental representation. As is well known (see, e.g., ||12||), various forms of the 
energy-momentum-tensor operator are available. Although they lead to different 
stress-energy-momentum densities, the same total or integrated stress energy mo- 
mentum results. Since we are interested in the stress-energy-momentum density of 
the quark-gluon plasma, it is of vital importance which form of the energy-momentum 
tensor is the physically sensible one. A natural candidate is the gravitational energy- 
momentum tensor M ^, since its generalization to the case of curved space-time ap- 



pears in the Einstein equation. As in all other known cases, Ml/ coincides |]T2[ with 
the Belinfante tensor: 

= e<g + e<g, (2) 

f\(g) — _P« pap 4. I n paaPpa 

+ Q a ){D v u) a + \A^d v d ■ A a ) + {// <-> v] 

+g^[{d p Q a ){D p uj) a 
A 
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(d ■ A a Y - \A a p d p d ■ A a , (3) 
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= \^ (jn Du +lu D^j i> - ig^A $ ip , (4) 



where 



(D^y^d^ + igf^Ay, 

Here {f abc } are the structure constants of su(N) and {T a } are a ant i- Hermit ian basis 
of the fundamental representation. The temperature-dependent part of the energy- 
momentum of a quark-gluon plasma is obtained through 

Tr (e~ H l T Q 



as 

p M = (e 0A .) - <o|e OM |o> , (6) 

where 

<o|e^|o> = lim(e^). 

It is to be noted that, in @, the ultraviolet (UV) divergence involved in (0o^) and that 
in (O|0o M |O) cancel out. Similarly the temperature-dependent part of the hydrostatic 
pressure P is defined as 

P = (e u )-(o|e u |o). (7) 

We employ the imaginary-time formalism [|l^] of thermal field theory. Never- 
theless, all the formulae are displayed in the Minkowski metric. Following 0, we 
introduce an index "r" , r = + for bosons and r = — for fermions: 

— = A r (P) = ( A+(P) ' f ° r b ° SOnS ' (8) 
P 2 I A _ (P) , for fermions. 

Here, for r = + { — }, po takes the values po = 2ninT {(2n + l)iirT} with n = 
0,±1,±2, ... . Capital letters represent four-momenta, lower-case letters their com- 
ponents: P M = (po, p) with p = pp, and P 2 = p$ — p 2 . For a loop momentum P, we 
introduce the symbol Tr defined as 

Tr F(p ,p)=T / jA- 3 F(Po,p)- (9) 

P n=-oc J ( 27r ) 
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For bosonic {fermionic} P, po = InmT {(2n + 1)ittT}. Tr is defined as in (0), 

soft P 

where the integration over p is carried out over the soft-p region, Eq. ([I]). 
In this paper, we use Feynman gauge (A = 1) throughout. 



3 Contributions from tadpole diagrams to (0^) 

To the lowest order, (|5|) is diagrammed in Fig. 1: Figs. 1(a), (b), and (c) represent 
the contributions from soft gluons, FP-ghosts, and quarks, respectively. Incidentally, 
for hard gluons [quarks], the HTL-resummed effective gluon [quark] propagator - 
the line with a blob — in Fig. 1(a) [(c)] can be replaced by the bare counterpart. 
The form of their contributions are well known, and we do not reproduce them here. 



3.1 Gluon sector 



We start with the gluon sector, Figs. 1(a) and (b). Using (^), we obtain for the 
contribution to (0^), 



Fig. 1(b) 



<©? ) = <eg? ) Pig iw + <eg?> 

1(a) = ~ (N 2 ~ 1) Tr p 7W(P) (#(P) , 



(©if) 



Fig. 1(b) 



N 2 - l) Tr A+(P) U V P 2 - 2P /1 P U 



soft P 



where *A Q/ g(P) is the effective gluon propagator in Feynman gauge, 

*A a ,(P) = E 

and, in (|TT|) , 



=T,L 



p 2 - MP) ' 



■P,P v g aP -P 2 5^ v + \g, v P 2 g^) 



In (HI), 



Q { S = E god gpj (Pi Pi - <y , 



(P = P/|Pl 



(10) 

(11) 
(12) 

(13) 



(14) 



5 



Q% — 9a0 9/30 — 



PaP/3 



— ^2 9 at 9/3jPi Pj , 



(15) 



and 11^ (n^) is the the transverse (longitudinal) part of the thermal vacuum po- 



larization tensor ||13|| , whose explicit forms are not necessary for the purpose of this 
paper. 



3.2 Quark sector 

Fig. 1 (c) with (Hj) gives 



-N f Tr 

2 soft p 



{Putr [ 7 /#(P) 
-g^tr [j?*$(P)]} + {ia++u} 



where *S(P) is the effective quark propagator [|14 



*S(P) = Y,^a*S (cT) (P) 

cr=± 

Pa = (l,ap). 



(16) 



(17) 



The explicit form of *S^ a \P) in ( |TTD is given in [TJ 



4 Hard- thermal- loop resummation for diagrams with 
operator 0^ insertion 

4.1 General observations 

Now we analyze formally 0(g 2 ) contributions to The relevant diagrams are 

depicted in Figs. 2 and 3. As stated in Sec. I, our analysis goes from the "diagram- 
matic point of view": We apply the HTL-resummation scheme |l], |2], |3|, |]] to each 
diagram in Figs. 2 and 3, and extract the contributions at high-T, which are of the 
same order as those from the lowest-order diagrams, Fig. 1. 

In each diagram in Figs. 2 and 3, the part enclosed with the box is the one-loop 
correction to the composite vertex, which we write (P|0 M „|P). We shall see that, for 
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hard K (~ T) and soft P (~ gT) [cf. Figs. 2 and 3], some of these matrix elements 
are HTL's, i.e., are of the same order as the lowest-order counterparts, and some of 
them lead to the same-order contributions to the energy and pressure densities as 
those from the lowest-order diagrams. 

The HTL's summarize the leading thermal effect and are essentially classical ||15|| . 
Then, in particular, they are free from UV divergences, so that we do not need 
to take the operator-mixing problem into consideration. In this relation, it is worth 
recalling here that the T-dependent parts of the above-mentioned (P^^P) represent 



|2|. [4|, |TT| essentially the matrix elements in tree approximation in vacuum theory. In 
order to see this in the present case, we take massless scalar </> 3 theory, for simplicity, 
and consider the composite operator <ft 2 /2. We compute one-loop correction to the 
composite vertex as depicted in Fig. 4: 

<Pl \ 2 (O) \V) = -y Tr {A + (K)} 2 A + (K + P) . (18) 

The manipulation of (|1^) is carried out in Appendix B. After continuing p (= 2ninT) 
to real energy p + i0 + , we have (cf. ( |B.4| )) 

(P\\<P 2 (0)\ P ) = \? j DuWD^K) 

xD n (K + P), (19) 

the formula which is written in terms of the real-time formulation of thermal field 
theory |12[, which is formulated on the time path — oo — > +oo — > — oo — > —oo—i/T, in 
a complex time plane. In (|T9|), P and K are the four- vectors in the Minkowski space, 



"1", "i", and "j" are the thermal indices, and Dij = 1,2) is the -component 
of the matrix propagator of a massless scalar particle, whose explicit form is given, 
e.g., in [|l|, [16[ . 

A thermal amplitude essentially represents the sum of the corresponding matrix 
elements J2ij(i\S\f)(f\S*\i), where S is the S'-matrix in vacuum theory, S* is the 
complex conjugate of the S'-matrix, and \i) [\f)] is the initial [final] state, which 
includes all the particles in the heat bath. The details and the rules for the corre- 
spondences are presented in We can apply the rules to each term in ([19]). As an 
example, we take ([19|) with i = 1 and j = 2 with k Q > and k + p > 0: 

1 



8tt 2 



K 2 + i0+ 



5(K 2 ) ®5((K + Pf) n B (k) (1 + n B (E)) 
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{2m 5 2 (K 2 ) <g> 8((K + P) 2 )} n|(fc)(l + 



(20) 



where £7 = |k + p|. It should be remembered |T2L ITBJ that, although products of 



singular functions appear in (0), they are cancelled in (|19D, leaving well-defined 
terms. The first and second terms in fl2~0|) represent the products of tree amplitudes, 
as depicted in Figs. 5(a) and (b), respectively. In Fig. 5, the left-side part of the 
final-state cut line (dot-dashed line) represents the amplitude or S'-matrix element 
in vacuum theory, while the right-side part represents the complex conjugate of the 
amplitude. The state at the left end and the state at the right end are the same and 
they represent the initial state. The lines on top of Fig. 5 stand for the spectator 
particles, which are the constituents of the heat bath. The left-side part of ® in 
the first [second] curly brackets in ( pOp stands for the amplitude S in Fig. 5(a) [(b)] 
and the right-side part stands for complex conjugate of the amplitude, S*, in Fig. 
5(a) [(b)]. Then, each term in (^) is simply an integration over k of the product of 
the tree amplitude and the complex conjugate of the tree amplitude weighted with 
statistical factors. Since the latter factors dump at the UV region, no UV disaster 
arises. This observation applies to all other portions of flTD|), provided that the T = 
part has been subtracted. 

The above argument applies to the T-dependent parts of all (PlQ^P)^ in Figs. 
2 and 3. Since no UV divergences appear, (P\Q 0fl \P) [(P\Q U \P)} enjoys a meaning 
of energy- momentum [pressure]. 

To summarize, as far as the T-dependent part of (6/^), which yields leading- 
order contributions at high-T, is concerned, we do not need to worry about the 
operator-mixing problem, and we can use the operators (|3|) and for evaluating the 
energy-momentum density. 

Now we separately study the soft-gluon case, Fig. 2, and the soft-quark case, Fig. 

3. 



4.2 Soft-gluon sector 

In this subsection we extract from formally 0(g 2 ) diagrams, Fig. 2, the contributions 
which are of the same order as those from the tadpole diagrams Figs. 1(a) and (b), 
Eqs. © - ©. 
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The contribution of Fig. 2(a) is written as 



<ew> 



Fig. 2(a) 



x*A 7<5 (P) ^(K)V^(i^, P) V e Sa (K + P, -P) 



(21) 



where 



V^(K, P) = g Py (P - K) t - g je {K + 2P)„ + g eP (2K + P) 7 . 
In (21])) *A 7 ,5(P) is as in fljjD and 6^(K) in (0). Carrying out the contractions with 



respect to the repeated indices and changing the integration variable K to —K — P at 
several places, (|2~T|) is led to a sum of terms, which are of the following generic forms, 



Tr — — P(P, K) <g> m 

softp p 2 -n^(P) v ' ! l] 

Tr -4: G(P, if) ® I++ , 

soft P P 2n y ' %3 



T,L) 



in 



1,2), 



(22) 
(23) 



where 



P(P, K) ® I, 



rr 



Tr F(P,K) {A r (K)Y 



K 



{/Y\K + P)} 3 



(r, r' = ±1 . 



(24) 



Keeping the terms which lead to the contributions of the order under consideration, 
we have 



<e<£> 

where 



g 2 N(N 2 - 1) Tr 

V ' soft p 



E 

!=T,L 



- y w + — v (4) + — v (2) 

p2 _ IL>(P) ^ P 4 ^ P 2 



+ (K\Q^\K) + cf (K^l + K V Q$) K> + cfK.K v 

+ (4V \k\QP\k) + cf (K^l + K V Q$) K p + cf K ,K u 
+ cfK,K v (K\Q^\K)®I+t, 
VjS = dxP^PJta + P-K {d 2 (K fl P u + K V P^) + d 3 (P ■ K)g^) ® /+ , 



(25) 



J 20 

(26) 
(27) 
(28) 
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In (gSp - ([§), Qfl is as in (|TJ) and 

(k\QM\K) =K p K a Q^, (29) 

cf^cf* for j = 2, 3, 4, 6, 7, 9 , (30) 

cf^cf) for J = 1,5,8. (31) 



The coefficients, c's and cfs, in ( |25|) - (|28|) are listed in the first row of Table I. 



In a similar manner, we can show that the contributions of Figs. 2(b) - (f) may 
also be written in the form (E3) - (EBp with (RDf) and fl3~T|). [The dashed-lines in Figs. 



2(e) and (f) stand for the hard FP-ghost propagators.] For each contribution, the 
coefficients are tabulated in Table I. 

Now we turn to analyze Figs. 2(g) - (i), where the dashed lines stand for the 
FP-ghost propagators. One ghost-gluon vertex is proportional to the soft ghost mo- 
mentum. Then, this ghost-gluon vertex brings in one power of a soft momentum, 
instead of a hard loop momentum as in the case of Figs. 2(a) - (f). Then, the 
contributions of Figs. 2(g) - (i) are negligible, 

(Q^)^0. (32) 

Finally, we analyze Figs. 2(j) and (k). It can easily be shown that, adding the 
contributions from Figs. 2(j) and (k), the u g fJiU term" in (U) - - the term that is 
proportional to g^ v — leads to vanishing contribution. It is worth pointing out in 
passing that the g^ u term vanishes provided the equation of motion is used. Then, 
what we should analyze is Figs. 2(j) and (k) with the first term on the r.h.s. of (^). 
It is an easy task to see that the term P a Pp/P A of the effective soft-gluon propagator 
*A a p(P) (cf. (|i~3D) does not yield leading contribution. Ignoring all other terms that 
lead to nonleading contributions, we obtain for the contribution to (6$), 

<^>U = ^(^-ISpl^idi^f* 8 ^ 

+c^K,K u ® IJ - 2 (jCfijg + K ,Ql2) K° ® I~ 
+4K lt K v (K\Q<- £ >\K)®In} , (33) 

a (T) =2a (L) = _ 2 _ 

In (|3"3"|), I^f e ^ c - are defined as in (El). 
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4.3 Soft-quark sector 

In this subsection we extract from formally 0(g 2 ) diagrams the contributions which 
are of the same order as that from the tadpole diagram Fig. 1 (c), Eq. (|16|). 

The diagrams to be analyzes are depicted in Figs. 3(a) - (c). The leading contri- 
bution from Figs. 3(a) and (b) is obtained as 



<e<g> 



quark 



2g 2 N f (N 2 -t 
x S S E (P) {4K»K U (K ■ P a ) ® I 2 + 

<T=± 

+ (p^ + P ffi/ i^) ® Ji + 



while Fig. 3(c) yields 



<e<g> 



gluon 



-±g 2 NAN 2 - 1) Tr V *S (ct) (^) 

soft P ' 



<T=± 



2K fi K v (K-P n 



K P + K P 



5 Energy density of soft modes 

In this section we compute the energy P , Eq. (|G|), of soft modes. 



(34) 



(35) 



5.1 Gluon sector 

Tadpole diagrams, Figs. 1(a) and (b). 

From ( |T0D - ([15]) , it is straightforward to evaluate the energy of soft modes, 



»G?o) 



(cf. (^|)). Employing the spectral representation ( |A.1| ) in Appendix A for 
(P 2 - n^(P)) -1 in (|T3]) and using flCO]) in Appendix C for A+(P) in (|T2j), we ob- 
tain 



<efe 0) > 



1 / \ r 1 l T r+°° 

2prK) (?5 + P 2 ) +PL«) (f 



, 2 2 



~ N 



1) Tr 1. 

soft P 
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where 



n B (0 



1 



e «/r - 1 ' 

and the spectral functions, p's, are defined as in ( |A.2| ). Summation over n (cf. 
and integration over r lead to 



(36) 



(6 



oo / 



N 2 - 1 
N 2 - 



d 3 p 
(2^)3 

Tr 1 

soft P 



P 2 )PL 



(£)} 

(37) 



Here and in the following, unless otherwise indicated, integration over p is carried 
out over the region (JT|) . In deriving fl3~7|), use has been made of the spectral sum rule 
( |A.3j ) in Appendix A. It is to be noted that, in Pq 9 °' ) (cf. (|6])), the last term in (|37j ) is 
cancelled by the term coming from (0|e& 0) |0). We recall that we are interested in the 
region p < O(gT) and note that Pt/l(0 vanishes |L3| in the region |£| > p + O(gT), 
where the term of O(gT) comes from the pole of {P 2 — Ui(P)) . Then, in the region 
of our interest, £/T << 1, and we may approximate ng(£) as 



T 

Using the sum rules, ( |A.3[ ) and ( A.4[ ) in Appendix A, we obtain 



,(90) 



(e^ 0) ) - 

(N 2 - 1) 



if 10) 
d 3 p 



(0|O 00 



(27T) 



2T + 

+ o 



m 2 T 



2 p 2 + 3 m? 



p 
T 



(3? 



(39) 



(40) 



where m s is the effective gluon mass induced by the quark-gluon plasma, 



ATA g 2 T 2 
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(41) 



In (|40|) , — p in the square brackets has come from — <0|O^ 0) |0> in ©. If we were to 
use the bare thermal gluon propagator g a/3 /P 2 instead of *A a/3 (P) in fllTD, we would 



obtain (|40"D with the second term in the square brackets being absent. It is to be 
noted in passing that (EH) is valid for T » p. Then, it is of no surprise that 
does not vanish in the naive limit T — ► 0. 
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Figure 2. 

Figs. 2(a) - (f): The starting formulae are - Q3"T| ) with Table I. All the 
necessary formulae for evaluating (6 o) are displayed in Appendix C. The result may 
be written in the form, 



;luon 



g 2 N(N 2 - 1) T J 



21 
16 



g 2 N{N 2 



1)T 



d 3 p 
(2^)3 

d 3 p 



7 T 2 



T 2 



16 p 2 16 p 2 + 3m 2 



m 2 g T 2 



(2ir) 3 p 2 (p 2 + 3m 2 



(42) 



which comes from small parts of the integration region in (PT|), where K is hard, 
k = 0(T) [HTL]. In (H), the contribution to the coefficients 7/16 (-7/16) from 
Figs. 2(a) - (f) are, in respective order, as follows, 



7 _ 49 11 1 3 3 1 

T6~48 - 24 + 8 _ T6~T6 + 8' 
7 1113 11 

_ 16 ~ ~ 16 ~8~8~16 + 48 + 24' 
Figs. 2(g) - (i): As discussed in Sec. IV, in conjunction with (|32|) , the contributions 
are nonleading as compared to (|40|), in which the first two terms in the curly brackets 
are kept. 

Figs. 2(j) - (k): The relevant formula is (^). Using ( |C.23| ) in Appendix C together 
with flCj) - ( EUD , we obtain 



quark 



\g 2 N f (N 2 



-IT 



T 2 



d 3 p 

(2tt) 3 p 2 + 3m 2 



(43) 



Adding all the contributions ([40"D, ([42]) , (|32|) , and (]43]) , we obtain for the leading 



contribution to P 



(9) 



) 



»(</) 



<eg?)-<o|eTO 



IT 



(2vr) : 



2 + liV ^ 
16 I p , 



1 



2 (p 2 + 3 



rri; 



3 m. 



2 V4 



- {-N + N f UgT) 
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d 3 p 

softp (27r) 3 



(n 2 - i) t 



63 



N 



3 13N + 8N 



i 



8 2N + Nf \ p 



8 2N + N f p 2 + 3m 2 n 



(44) 
(45) 



where we have used (ffl]) . It is an elementary task to show that £ s (p) is positive 
definite for any values of N, Nf, T, and p. 

5.2 Quark sector 

Tadpole diagram, Fig. 1(c). 

The starting formulae are flTE|) with /i = z/ = and (]T7|) . Using the spectral 
representation ( A.5 ) in Appendix A, summing over n (cf. @), and carrying out the 
integration over r, we have 

d 3 p 



<eif> 



(2^ 



■p J d£{l-n F (0} 
E ^(0 



.<T=± 



where 



MO 



1 



(46) 



(47) 



e«/ T + 1 ' 

and the spectral function p a is defined as in ( |A.6| ). Since p a {0 is nonzero only in the 
region |£| < p + O(gT), we may approximate n F (C) as 



np(f ) — — + .... 



(4£ 



Then, using the sum rules ( |A.7|) - ( |A.10| ) in Appendix A, we obtain 



(©if) = # 



<i 3 p p 2 



where 



m. 



soft P (27r) 3 T 
1 iV 2 — 1 2 



1 + 



> 2 + 5m 2 /3 N 



(49) 



(50) 



J 16 

Note that the effect of the HTL resummation for the quark propagator appears at the 
term of 0(m 2 /T 2 ) = 0(g 2 ) relative to 1. In fact, with S(P) in place of *S(P) 
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leads to flUf) in which p a (£) = §#(£ — &p): 



(©if) 



-2JV, 



AT* 



ci 3 p p 2 
(2tt) 3 T 



Now we rewrite (|5~H) as 
(51) = 



4iV f 



(2tt) 3 
(2^ 



pn F (p) + (0|9if|0) 



Using (H), d|), and (||), we obtain 



»(«o) 



(©oo 



soft p (27r)' 



oo |0> 



£ q (p) 



N f 



2p + 



■£q(P) > 

V + 5mJ/3 N 
T 



(51) 
(52) 

(53) 
(54) 

(55) 
(56) 

(57) 
is 



In contrast to the case of gluon sector, the T-dependent part of £ q (p), Eq. 
nonleading. The origin of this difference may be traced back to the different statistics, 
(BBT) and 



Here some comments are in order. The factor "4" in (^) - ( |54| ) represents four 
degrees of freedom, 4 = 2x2, where the first "2" comes from a quark and an 
antiquark and the second "2" from spin degrees of freedom. pnp(p) in fl53|) represents 
the thermal energy per one mode of momentum p and — p/2 is the zero-point energy. 
In fact the zero-point energy per one fermionic mode of momentum p is obtained as 
(0|| (aj,a p — a p a^ p )\0) = — |, where a p [aj,] is the annihilation [creation] operator of 
an (anti) quark with p and with definite spin. It is interesting to note that, for soft 
p, np(p) ~ 1/2 (cf. (f|8[)) and then, in (|53"D , (the leading part of) the thermal energy 
and the zero-point energy cancel out to yield fl5*2"|). S q (p) in (|5"T| ) is the thermal energy 
of (anti)quarks with momentum p, Nf x 4 x p x uf{p) — 2Nf x p. 
Figure 3. 

Figs. 3(a) and (b): The contribution to the energy density is obtained from fl34| 



as 
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quark 



Ag 2 N f (N 2 - \ j 

xTr ^ *S (a \P) [2k 2 (K ■ P a ) ® 1 + (3A; - 2(xk ■ p) 



cr=± 



(5f 



With the help of (|A.5|) in Appendix A and ( |C.1| ) in Appendix C, it is straightforward 
to show that the 1 11 part in ([58]) leads to a nonleading contribution. As to the 
contribution with P ar t i n (f>B"D, using ( CTT ) and ( |C2| ), the relevant HTL parts 
are easily evaluated: 



T *S {a) (P)k 3 F(k-p)(g)I 2 l 

n=— 00 

1 /"+00 



1 f +Q ° ly- fs-S f d 3 k ^.f A 



X 



E 



=± 



— ^f)(1 + 2ub 



+ n B (l + n B )(l-2n F ] 



k • P - s£ 
sk • p 



(k- P -*0 2 J 



1287T 2 
/■+00 



C(3)T 2 | d(k-p)F(k-p) 



> 



+ 



sk ■ p 



+00 

T *S^(P)k k 2 G(k-i>)®I 2 + 

n=—oo 

7 



1287T 2 

r+00 



C(3)T 2 | d(k-p)G(k-p) 



x 



/+00 
-00 1 



+ 



k P 



±l_k-p-s£ (k-p-sO 2 . 



We note that, in (58), 



if ■ P a = k — ak ■ p 



(59) 



(60) 



(61) 



Eq. ( |58"D with the k part in ( pT| ) leads to ( |6U|) with G(k • p) = 1, which vanishes 
because the integrand is odd with respect to k • p. The other part in (|6lD , — ak ■ p, 
leads to ( |59"D with F(k ■ p) = — ak ■ p, which also vanishes due to the same reason as 
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above. Hence the contribution fl58|) is of nonleading 



<e$> 



, ~o. 

quark 



Fig. 3(c): The starting formula is ( 35Q . The same analysis as above shows that 



~ o. 

gluon 



In conclusion, to leading order at high-T, P^°\ Eq. Q5B|), does not receive addi- 
tional contributions from Figs. 3(a) - (c). 

6 Hydrostatic pressure 

In this section, we compute the hydrostatic or kinetic pressure density of soft modes. 
The T-dependent part of the hydrostatic pressure P is defined as in (^): 

P = (0 11 )_(O|0 11 |O) 

= -^[{(0, M )-(o|e/|o)}-{(e oo )-(o|e O o|o)} . 

(62) 

We are interested in the leading contribution at high T to (|62"D for soft modes. Then, 
as in Sec. V, no diverging integral appears. 
We shall show that, to leading order, 

(e/)-(o|e/|o)~o. (63) 

Using this in (|62"D , we obtain 

p^i{(e o)-(o|e o|o)} = ^Po, 

where Po is as in (HD with /i = 0. Then, defining the "pressure density" V(p) as 

Jsoftp (ZTT) 6 

we have, for the gluon [quark] sector, V g (p) — S g (p)/3 \P q {p) — £ q (p)/3 ], where S g 
[S q ], Eq. ( f£5|) [Eq. (|57|)1, is the energy of soft gluon [quark] with momentum p. 
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We are now in a position to confirm (|63|) . From with A = 1 and (|j), we obtain 



eto)" = 2 [(<9 M u a ) (D» uf - d» (Al d ■ A a 



0(ff)A» = -3z^ 



(64) 
(65) 



We note that a term with a total derivative does not contribute to (|63|) . Then, the 
second term in the square brackets in (|64]) may be ignored, and we may replace ( |64"D 
by 

= — 2cu a (9 A1 (D fM u) a . (66) 

It is to be noted in passing that (|65D and (|66| ) vanish when the equations of motion 
are imposed. 

The diagrams to be analyzed are Figs. 1(b), Figs. 2(e), (f), (g), and (h), with 
q(s)a* insertion, and Fig. 1(c), Figs. 2(j), (k), Figs. 3(a), and (b), with insertion. 
Fig. 1(b): We obtain for the contribution to (0^ M ), 



Fig. 1(b) 



-(N 2 -l)Ti^-- P 2 

V > p p2 

-(N 2 - l)Tr 1. 



(67) 



As in (|7|) above, (|7]) is cancelled by (0|9^|0) in (H): 



[(© { f ) - (o|e«"|o> 

Fig. 1(c): The contribution to (O^) is 



Fig. 1(b) 



0. 



(e ( f) 



3A^/ Tr tr 



*#(P) 



(6? 



Fig- 1(c) J soft P 

Computation of ( |68| ) goes through the similar procedure as in Sec. VB. In place of 
(BED, we have 



(6(f) = 12N,Jdt{l-n F {Z)} 

d 3 p 



Then, using (^8|) and ( |A..7| ) - (|A.10|) in Appendix A, we obtain 

d 3 p 



[(9(f) - (o|e(f |o> 



Fig. i(k) j soitp (2tt) 3 



o { t (p) 



( ? (p) ^ 



T 



(69) 
(70) 



where m/ is as in (|50|). 

In the region of our interest, Eq. ([I]), ( p9"| ) - (|7^) are subleading as compared to fl5"7|). 

Figs. 2(e) and (f): It is straightforward to show that the contribution of Fig. 2(e) 
cancels the contribution of Fig. 2(f). 

Figs. 2(g) and (h): The contribution is 

-2, 2 iV(iV 2 -l)Tr p ^/^. (71) 

Comparing ( [71] ) with the u cIq term" in fl28|) and noting that K (P) is hard (soft), we 
see that ( |7T|) is subleading. 

Figs. 2(j) and (k): It can easily be shown that the sum of the contributions of 
Fig. 2(j) and Fig. 2(k) vanishes. 

Figs. 3(a) and (b): The contribution is 

12^ 2 N F (iV 2 ~ 1) Tr £ £M(P) (A* - ak • p) I++ . 

soit P — 

<T=± 

This is of the same type as ( f)8D with the second term in the square brackets and is 
nonleading as compared to 



This completes the proof of fl63|). 



7 Summary and discussion 

Applying the HTL resummation scheme to the one-loop correction to the composite 
vertex (P\Qfj, u \P), we have deduced the energy density and the pressure density, in the 
high-T limit, for soft gluons and soft quarks in a hot quark- gluon plasma. The resul- 
tant energy density for soft gluons, £ g (p), and for soft quarks, £ q (p), are given in ( K5[ ) 
and (|57D , respectively. In contrast to the soft-gluon case, to leading order at high-T, 
the energy density for soft quarks, £ q (p), does not receive contributions from HTL's. 
This difference originates in different statistics. To leading order, the hydrostatic 
pressure density V(p) is related to the energy density £{p) through V(p) — £(p)/3, 
for both soft-gluon and soft-quark sectors. Thus the relation that is characteristic of 
free-massless-particle systems still holds approximately. 

Now, as in Sec. IV, we discuss the physical content of (@oo) by decomposing it 
into the sum of matrix elements of S <8> S*. We take up ( |TID and (pip, and expands 
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*A aP (P), Eq. (0), in powers of IL e (P): 



A afj {P) = *Ag(P) + E*A$(P) 



* A ( n ) / 



(72) 



n=0 



*A%\P) = P a P?{A+(P)} 2 , 



Inserting ( [72] ) into ([□]) and (PJ), we have, with obvious notations, 



, n+l 



(e? 0) ) 



and 



(©oo) 



(©?o 0) )F+E(©r) 



Fig. 1(a) 



Fig. 2(a) 



(90) \ 



n=0 



(©00)F + E( e 00)n 



(73) 



(74) 



n=0 



respectively. According to the general arguments in |18| (see also [19), through 
analytic continuations to the Minkowski space, ( |73|) and ( |74|) go to the expressions 
written in the real-time formulation: 



(eSfV = -\ (iv 2 - i) / ^ (P; a 2 ) 

(©?0 0) )n = 4 " !) / C(P) E 2 



A = 



(75) 



(0 

a/3 



J2(-y- 1 D l3 (P){[tle(P)r 3 D(P)} n } n , 
d A P r d 4 K 



(76) 



(©oo)f = -|iV (iV 2 - l) J — / _ C WV, 7 e(^ P) V fa (# + P, -P) 



x E D tl (K)D lj (K)D Jt (K + P)P,P s -^D lJ (P ] X 2 



A = 



(77) 



(e 00 ), 



' 2 



N ( N2 ~ 1 )J^Jj^ e^(K)v^(K, p) v fa (* + p, -p) 



(2tt) 4 7 (2tt) 



x E DhWDuWD^K + P) E c$5 

x E(-K- 1 p^(p){[n,(p)r 3 p(P)] n }^ . 

j'=l 



(7* 
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Here, as in (0), P and K are the four- vectors in the Minkowski space, D is the scalar 
2x2 matrix propagator, Tle(P) is the scalar 2x2 vacuum-polarization matrix, and 
r 3 is the third Pauli matrix. 



Applying the general rules [17[ to each term in (|75|) - (|7q), we can identify the 
corresponding set of matrix elements of S ® S*. As an example, we take up (Oqq ^ 
in ( |7B|) and (600)3 i n ([^), and find the corresponding matrix elements. Thermal 
amplitudes representing (Boq ^ an d (©00)3 are depicted, respectively, in Figs. 6(a) 
and (b). In Fig. 6, P is soft, K's are hard, and "1", "ii" — "V' are the thermal 
indices. We consider the case where i\ = 12 = 23 = i§ = 1, £4 = 25 = 2, and all k Q 's 
are positive. The rules ]T7|] tell us that the thermal amplitude, Fig. 6(a) [(b)], under 
consideration represents S <S> S* as depicted in Fig. 7(a) [(b)]. 

In S (g> S* in Fig. 7(a), the "probe" Boo is put on the soft line P in the S- 
matrix element side, while in Fig. 7(b), Boo is on the hard line K. From the very 
definition (|5|), (Boo) is the (thermal) expectation value of Boo the initial states - 
the states at the left end and the right end of Fig. 7, which consist of hard active 
gluons only. Here "active gluons" means those gluons which directly participate 
in the thermal reaction. We recall here that, as has been stressed by Weldon p0[ 
in a different context, intermediate states or virtual particles are a mathematical 
figment of Feynman-Dyson perturbation theory. Then, Fig. 7(a) [(b)] represents the 
contribution through the figment, intermediate soft [hard\ gluon with P [K], to the 
energy density of the quark-gluon plasma, which is composed of interacting gluons 
and quarks and is in thermal equilibrium. Thus Fig. 7(a) [(b)] has nothing to do 
with the energy of the soft [hard\ gluon. As a matter of course, Fig. 6 or the original 
Figs. 1, 2, and 3 include the diagrams, S <S> S*'s, in which Boo is put on the gluon 
line in the initial state. 

In this relation, we may view Figs. 2 and 3 as giving the correction to the energy- 
momentum density of hard gluons and quarks. To the lowest order, the energy density 
of hard gluons is obtained from (|3"7D by setting pr(Q = Pl(£) = <K£ 2 — P 2 )i 

P ( o 90) =2(iV 2 -l) f£!Lkn B (k). (79) 

hard gluon J [ZTT) 

The contributions from Figs. 2(a) - (f) to the energy density are of the following two 
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types: 



2 „ f d 3 p 1 f d 3 k 1 
9 T J tfi + (3m;) 7 * W ' (80) 



(2tt) 3 p 2 + (3m 2 ,) y (2tt) 3 jfe 
(2tt) :{ /r' 4 (3m 2 ) 7 (2-) : < 



^1) 



Carrying out the integration over k in each contribution, and summing over the 
contributions from all other diagrams, we obtain (|44|) and (f45|). In order to obtain 
the correction to P^ , Eq. (|79|), from Figs. 2(a) - (f), we carry out the 

hard gluon 

integrations over soft-p in (|80| ) and (pTj) : 

Eq. (80) oc g* p *T J -^±n B {k), (82) 

Eq. (H) oc 0Vy-^ nB (fc){l + n fl (A:)}. ( 83 ) 

Here p* denotes the boundary between the soft- and hard-P regions, gT << p* << T. 
Eqs. ([82] ) and ( |33"D is 0(p*/T) smaller than (f79|) in the hard k {= 0(T)} region, which 
gives the leading contributions to the integrals in (f79|), flS2]), and (|83|). Thus Figs. 
2(a) - (f) yield 0(p*/T) smaller corrections to the hard-gluon energy density. This 
is also the case for Figs, 2(g) - (k). The diagrams Figs. 3(a) - (c) yield higher-order 



corrections to P^ 



hard quark 



Appendix A Effective propagators 

In this Appendix, we display the spectral representations for the effective propagators. 
Effective gluon propagator (cf. (|13[)) 



1 

p 2 -n T/L (P) 

= - f^dTe^ [ +Q °dtp T/L (0(l + n B (0)e-^, 

JO J-oo 
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Pt/l{£) 



1 

-— Im- 



Po 



(A.l) 
(A.2) 



vr p 2 - n T/i (P) 

where %(0 is as i n (p6|)- Explicit expression of TIt/l(P) is given in [|n||. The spectral 
function, Pt/z(6)> defined above satisfy various sum rules, among which we use the 
following ones in the text; 



oo 
' oo 



d^ n PT/L(0=n, (n = 0,l) 

... 1 



f S. (n — 



(£ = T, L) 



In ( |A.4| ), m 5 is as in 



(A.3) 
(A.4) 



Effective quark propagator (cf. (|T 



*S^(P) = - / dre P0T / dCe^ r p CT (0 

x (1 - n F (0) (a = ±) 



M0 = — /m^)(P) 

7T 



PO = £ + i0+ 



(A.5) 
(A.6) 



where np(£) is as in (f47|). Explicit expression is given in (14| . The spectral 

function, p a (£,), defined above satisfy the sum rules, 



/+oo I 
-oo Z 

/+oo 1 
-oo Z 

/+oo I , 

-oo Z 



(A.7) 
(A.8) 
(A.9) 
(A.10) 



where m/ is as in 
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Appendix B Derivation of (19) 

In this Appendix, we deduce (0) from ([181) . Using the "r-representations" ( |C.lD and 
( |C.2|) (in Appendix C) for A + (K + P) and {A + (K)} 2 , after standard manipulations, 
we obtain 

.2 



- 9 -Tr{A+(K)YA+(K + P) 



16 



■°° d 3 p 1 



[1 + n] 



n 



l + cr 
2 

1 - a 



n 



n 



1 



x 



+ 



+ 



1 



k k k + oE + rpo ' (k + aE + rp ) 2 
nil +n) a 



T k + aE + rp 

where n = n^(fc) and n' = tib(E) with E = |k + p|. We note that 
written as 

g 2 d r d 3 k 1 



(B.l) 
may be 



8 d\ 2 J (2tt) 3 EE x 

1 + n B (E x )) (i±2 + n') - tib{Ex) (if + n> 



x E E 

T=± <X=± 



A = 



(B.2) 



where ^ = V& 2 + A 2 . In the form ( B.l|) or ( B.2|) , we can continue po (= 2mnT) to 
real energy po + z0 + . Straightforward calculation shows that 



(B.2) 



2 9 



diK h 

i=i 



(2vr) 4 



X 



^DaiK + P) 
dD u (K; A 2 ) 



<9A 2 



A = 



(B.3) 



The formula (P-3| ) is written in terms of real-time thermal 3 theory |12[ , which is 
formulated on the time path — oo — > +oo — > — oo — > — oo — i/T, in a complex time 
plane. Du(K; A 2 ) is the bare thermal propagator of a boson with mass A. In passing, 
the result (B.3|) is in accord with the general argument of analytic continuation of 
thermal amplitudes (cf. [[Lp, [H 
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Employing the mass derivative formula |TL], [L6| for dDu/dX 2 , we obtain 

(B-3) = y 2 I iKZU E H i+j Dii(K) 

x Di^Dj^K + P) . 



1 2 r d*K_ 2 
2^ 7 (2tt)* 



(B.4) 



Appendix C List of useful formulae 

In this Appendix we collect various formulae used in the text. For the purpose of 
deriving those formulae, we introduce the "r-representations" |^|, [ISj for the bare 
propagators. As in Sec. I, following ||, we use an index r, r = + for bosons and 
r = — for fermions: 



n 



n + (k) = ns{k) , for bosons 
n~{k) = np(k) , for fermions 



and introduce 



f + (k) = 1 + rn r (k) , /I = rn r (k) , 
g r (k) =rn r {k){l + rn r (k)) . 

Then, the boson and fermion propagators (cf. (|])) may be written as 

1 rV T 



Ar{K) = ~2k I dTek ° T §. ^(*) e "* fcr ' ^ 
{A r (K)Y 



'0 s=± 
,2 1 1 



2k 2 K 2 

rl/T 



+ ^r gr(k) l! drekor \- e ~ kT + ek " 

+ —j o Tdre^^sf s (k)e-^. 



(C.2) 



Formulae for evaluating Figs. 2(a) - (f) 
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As an illustration, we sketch how to evaluate, under the HTL approximation, 

+00 7,2 +00 1 

, , x (8> Itt = T Y ———— Tr k 2 A + (K)A + (K + P), 

(£ = T, L) . 



T E 



Using ( |C.1| ) for A + 's, standard manipulations lead to 



Tr k 2 A + (K)A + (K + P) 
1 r d 3 k k ^ 



1 + n B (k) + n B (E) 
k + E + rpo 
n B (E) - n B (k) 



k - E + rp 

where i£ = |k + p|. We take out from ( |U.3| ) the UV-divergent piece; 



(C.3) 



1 r d 3 k k 



^E 



4 J (2tt) 3 £ T ti fc + E + t Po 
Continuing the form ( |C.4j ) to real po, we can easily show that 

"+°° dk n r d 3 k k 2 1 



(C.4) 



—1 



-00 2tt J (2tt) 3 K 2 (K + P) 2 ' 

which is nothing but the "T = sector" of ( |C3| ) . Thus ( |C3[ ) is not the HTL and we 
ignore it. 

Using (|A.l ) for (P 2 — n^(P)) -1 , under the HTL approximation, we have 

+00 L2 



t E 



p 2 - MP) 

\ r+oa 
4 7-oo 



A? ^ — ^ 



k 



n B (k) n B (0 



-T- 



(1 + n B {k))n B {k) - t (n B {k) - n B {E))n B {£) 



p ■ k - r£ 

Since £, p << k, we can employ the approximations, 

n B (k) - n B (E) ~ n B (k)(l + n B (k)) 
T 

"b(0 - y , 



(C.5) 



p_k 

T 
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to get 



(C.5) 



1 [+°°dZ f 

2 L Y J 



d 3 k 
(2^)3 



n B (k) + n B (k)(l + n B (k)) 



Using the spectral sum rule ( A.4 ), we finally obtain 



T E 



k 2 



T 3 



p 2 - MP) 



%p 2 + ?>5 lL m 2 g ' 

(£ = T,L) , 



where m g is as in fl4"ip. 

Evaluation of leading-order contributions of other formulae used in the text are 
obtained in a similar manner: 



+oo 



T E 



/ w _L 

Jin — ~ 



T 3 



P 2 -U e (P)^ 10 12p 2 + 3^ L m 2 ' 



(£ = T, L) 



+oo 

T E 

n=— oo 



k 2 



P 2 - IL e (P) 

+oo r,2 
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r E 



p 2 - MP) 



i 



T 3 



1 1 



p 2 + 3 5 €L m 2 



+oo 

r E 

n=— oo 

+oo 

t E 

n=— oo 
+oo 

t E 



A:' 1 



p 2 - MP) 



01 



T 3 
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32 p 2 + 3 5, L m 2 



P 



P 2 -Tl L (P) P 2 



1 10 — 



T 3 



12 p 2 + 3 m 2 



A- 2 



p 2 -n L (P) p 2 



20 



+oo 



t E 



+oo 



t E 



p 2 - n L (p) p 2 

i (p-ir) 2 
p 2 -n L (p) p 2 5 



p i 

++ i 

1 ii 



T 3 



8 p 2 + 3 m 2 



1 T 3 

J 1 1 _ 



24 p 2 + 3 m 2 



+oo 

T E 

n=— oo 



1 (P-K) 2 , 1 T 3 

TT / D\ D2 2 r) /l „0 , o 



p 2 -n L (P) p 2 



24 p 2 + 3 m 2 



(C.6) 



(C.7) 
(C.8) 

(C.9) 



(CIO) 



(C.11) 
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±5? k-p PK ++ 1 T 3 



P 2 -n L (P) P 2 11 24p 2 + 3m 2 ' 



k-p P-if + 1 T 3 



n £- oc P 2 -n i (P) P 2 W 20 24p 2 + 3m 2 



m +^ k 2 (P-K) 2 5 T 3 

n £^ 00 P 2 -n L (P) P 2 ^ 21 96j9 2 + 3m 2 



deriving - ( EH) , we have used (jCJ) for 1/P 2 = A+(P). 



Using the above formulae, we obtain 



1 1 



«• [K\QF>\K) ^ 

L p2_n r( p) ®^o- 12p2 



n=— oo 



n=— oo 



2- p2-n T (p) ® 21 - 48 p 2 



P 2 -n L (P)^^ 11 24p 2 + 3m 2 



+oo 



+ 
20 



+^ feo^Qg + 1 T 3 



P 2 -n L (P)^ i20 24 p 2 + 3m 2 



P(K|Q( L )|ir) i t 3 
— ^ — _ g> Itt ~ — — - 



T JL P 2 -n L (P) ®^ 21 32p 2 + 3m 2 ' 



+oo 
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Here, as in @g), [k\Q^ t / l ^\K) = K" K° QfJ L) . 

For evaluating the formula of the type (|2"3|), further manipulations are not neces- 
sary. In fact, with the help of the trivial identity, 

+00 1 



p2 

n=— 00 



lim T 



xG(P,K)®If+, (£ = 0,1), (C.22) 
we can evaluate the l.h.s. of (|C22j ) by using (|C.6|) - ( |C.21|) . 



Formulae for evaluating Figs. 2(j) and (k) 

When I~~ is substituted for Ifj in each of the formulae above, (|C.6|) - ( |C.21|) , 
we obtain the r.h.s. of that formula multiplied by the factor —1/2. Namely, writing 
the l.h.s. of - ( P^lD in a generic form T£i~ H(P, K) <g> /++, we have 

+00 1 +00 

T J2 H(P,K)®l;r^--T J2 H(P,K)®lt+. (C.23) 
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FIG. 1. Tadpole diagrams. The symbol (g> indicates the point where the operator 
Q^ v is inserted. The blobs indicate the (HTL-resummed) effective propagators. 

FIG. 2. Formally 0(g 2 ) diagrams that contribute to the soft-gluon sector. K is hard 
(~ T) and P is soft (< O(gT)). 

FIG. 3. Formally 0(g 2 ) diagrams that contribute to the soft-quark sector. K is 
hard (~ T) and P is soft (< O(gT)). 

FIG. 4. One-loop correction to the composite vertex (P\<ft 2 /2\P). 

FIG. 5. A contribution to S®S*, which is included in the thermal Green function Fig. 
4. The left-side part of the final-state cut line (dot-dashed line) is the 5*-matrix 
element in vacuum theory and the right-side part is the complex conjugate of 
the S'-matrix element, S*. A group of particles on top of the figure stands for 
the constituents of the heat bath. 

FIG. 6. Diagrams representing the thermal amplitude (Ooo°^)3, Eq. (|76|), and (Ooo)3> 
Eq. (|78|). For visual clarity, we have used solid lines for gluons. Thick lines 
represent hard gluons, while thin lines represents soft gluons. 

FIG. 7. Matrix elements of S <S> S*, which are included in the thermal diagrams Fig. 
6. 
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Values of the coefficients in (26) - (28) for various diagrams in Fig. 2. 
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